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1. Introduction 

The phase diagram of strongly interacting matter shows an astonishing richness (cf. Qj], ||] 
for some excellent overviews). At finite temperature T CP ~ ^(160MeV) and non-zero baryo- 
chemical potential jj,jj p ~ <^(360MeV) Q3|, |] a first-order phase transition line (for fX B > /J.g P ) 
separating confined from deconfined matter terminates. While the exact location (T CP ,jUg P ) of 
this end point, the critical point (CP) being of second order, is matter of current investigations (see 
comparison in showing in fact a strong quark mass m q dependence ^]), its existence for 
finite quark masses seems to be guaranteed by general arguments [||, |8]] relying on universality 
arguments. Note, however, the recent discussion [§, [RJ, where an alternative scenario without CP 
is discussed. We take here the attitude to assume, in line with the existence of a CP and discuss 
phenomenologically some of its possible implications. Besides, the extension of the critical region, 
which might be crucial for searches of CP in future heavy-ion experiments, is fairly unknown 



[ ]11| , |12fl . For massive quarks, the CP belongs to the universality class of the 3-dimensional Ising 
model which was numerically verified by first-principle QCD simulations (dubbed lattice QCD) for 
the volume dependence of the Binder cumulant [EL O). In the (unphysical) limit of vanishing quark 



masses, the CP turns into a tricritical point belonging to the universality class of the Heisenberg 



model in three dimensions [ |14] ] while the transition line continues as second-order phase 
transition line of critical points up to \i% = 0. 

For finite m q and ji B < jj,g P , the sequence of first-order phase transitions turns into a cross over. 
Still, in this region a rapid increase of the entropy density s or energy density e in a narrow temper- 
ature interval or a peak in a suitable susceptibility allows for the determination of a pseudo-critical 
line defining the deconfinement temperature T c at Hb = 0. At T c , the chiral condensate < qq > as 
order parameter of chiral symmetry breaking and the expectation value of the Polyakov-loop as 
order parameter of deconfinement also rapidly change. For realistic quark masses reproducing the 
physical hadron spectrum, the transition of confined (hadron) matter to deconfined (quark-gluon) 
matter seems to be indeed a cross over [Q, |TJ]. This is in contrast to pure gauge theory, which 
displays at T c a first-order phase transition [[J, |i6fl . The quark mass dependence of the transition 
characteristics at T c is discussed in some detail in [|2, 15]. 

The CP as fundamental issue of QCD is matter of current investigations both theoretically 
and experimentally. In particular, observable consequences are discussed [i7[ |. Concerning the 
hydrodynamical description of the expansion stage in heavy-ion collisions, one might ask to what 
extent the equation of state (EoS) becomes modified by CP. 

Our paper is organized as follows. In section 2, we briefly review our quasi-particle model for 
the EoS of strongly interacting matter and compare with recent lattice QCD simulations for two 
quark flavors, in particular the isentropic evolutionary paths during the hydrodynamical expansion 
in heavy-ion collisions. In section 3, we discuss the influence of including critical point effects 
phenomenologically into our model on the EoS. Our results are summarized in section 4. 

2. Quasi-Particle Model and EoS for N f = 2 

The pressure p of Nf = 2 light quark flavors in our quasi-particle model (QPM) as function of 
temperature T and one chemical potential ji q = jUg/3 reads p(T,jj, q ) = Y, a =q, g Pa — B(H q ^ g (T , jj, q )) 
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(cf. [18] for details) with partial pressures of quarks, p q , and transverse gluons, p g , predominantly 



propagating on-shell with dispersion relation co a = \Jk 2 + Tl a . This quasi-particle picture is mo- 
tivated by a chain of feasible approximations starting from the Luttinger-Ward approach to QCD 



(cf. |20p for details). Yl a are given by the asymptotic expressions of the gauge-independent 
hard-thermal loop / hard-dense loop self-energies [ pi] ] and B is determined by thermodynamic self- 
consistency and stationarity of the thermodynamic potential under functional variation with respect 
to the self-energies, 8p/8Yl a = [|22|]. Non-perturbative effects are accommodated in the model by 
replacing the running coupling g 2 in II fl by an effective coupling G 2 which needs to be determined 
for arbitrary T and \i q . This can be achieved by solving a flow equation following from Maxwell's 
relation [|T8[ EO] as Cauchy-problem together with a convenient parametrization of G 2 at \i q = 



[23] 



From the pressure p, other thermodynamic quantities such as entropy density s or baryon 
density ns follow straightforwardly. Assuming local entropy and baryon number conservation 
during the hydrodynamical expansion in heavy-ion collisions, the evolutionary paths of individual 
fluid elements are represented by isentropic trajectories s/ns = const in the T - pLs plane. Recently 
riB, s or energy density e have been evaluated in first-principle numerical simulations of QCD 



[ |24| , |25[ ] as Taylor series expansions 

s{T^ B ) = T^s n (T)(^X , (2.1) 
e{T^ B ) = T^e n {T)(^X , (2.2) 



including terms up to order (hb/T) 6 (terms of order (/ig/T) 8 are currently calculated). Here, 
e n (T) = 3c n {T) +c' n {T) and s n {T) = (4 -n)c n (T) +c* n (T) where c' n {T) = Tdc n (T)/dT and 



c n (T)- l3 '^ T4 



(2.3) 

n q =o 



n\d(ll q /TY 

are the expansion coefficients for the pressure p [El] which can be derived from the QPM ex 



pression for p via (2.3). Using these expansions up to order (ju^/r) 6 , we calculate the isentropic 



trajectories for s/ns = 300, 45 and compare with lattice QCD results [|^] in Fig. 1. The agreement 



is fairly good (cf. [ [2q ] for a discussion of these results). In Q23|], we reported an impressive agree- 
ment between QPM and lattice QCD results for c n (T). As the Taylor expansion coefficients of s 
and e contain both, c n and c' n , they serve for a more sensitive test of our model. We compare the 



QPM results for ^2,4 and <?2,4 with lattice QCD results [25] in Fig. 2 finding a fairly good agreement. 



It should be emphasized that the QPM does not include any critical behavior. However, the 



adjustment to lattice QCD data for p(T,jj, q = 0) [27] forces us to build in a change of the slope of 
the effective coupling G 2 {T,ji q = 0): Above T c , the coupling resembles a regularized perturbative 
logarithmic behavior, while below T c a linear dependence on T is required to describe the data. 
This change of the slope of G 2 (T,/J. q = 0) at T c generates the structures seen in Fig. 2, as higher 
order derivatives of G 2 enter. With the above reasoning we now attempt to supplement our QPM 
by structures causing a critical behavior. 
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Figure 1: Isentropic trajectories: Comparison of Nf = 2 lattice QCD results [||| for s/ng = 300 (triangles) 
and 45 (circles) with the corresponding QPM results. In both cases, a truncated Taylor series expansion up 
to order (/Xb/T) 6 for s and «g is applied. The deconfinement temperature at jis — is set to T c = 175 MeV 

& 




Figure 2: Comparison of our QPM with lattice QCD results p5| ] for e n (T) (squares) and s n (T) (circles) as 
function of T /T c for Nf = 2; left (right) panel for n — 2 (4). QPM parameters chosen according to ^2% Efl ], 
The pronounced structures in the vicinity of T = T c follow from a curvature change in the parametrization 
ofG 2 (r,Ai,-0) (cf. ||). 



3. Including the Critical Point 



Following Q28p , we aim at incorporating phenomenologically CP features into the EoS 
parametrization and choose as starting point the decomposition of the entropy density s into a 
regular part s reg and a singular part s s j llg which is related to phase transitions and critical phenom- 
ena [|^]. For s reg we use our QPM entropy density, and s S i ng is constructed from the parametric 
representation [28, [M]] of the Gibbs' free energy density G(r,h) with critical behavior belonging to 



the universality class of the 3-dimensional Ising model. The order parameter of the 3-dimensional 
Ising model is the magnetization M(r,h). G(r,h) and M(r,h) are functions of reduced temperature 
r = (J — T c ) /T c and external magnetic field h. From p8l 150], the parametric form of G in terms of 



4 



Quasi-particle perspective on critical end-point 



Burkhard Kampfer 



new parameters R and 6 reads 

G = hoMoR 2 ~ a g(d)-Mh (3.1) 

with M = MqRP G and Mo, ho are normalization constants. The variables of the 3-dimensional Ising 
model are defined by 

r = R(l-d 2 ), (3.2) 

2 

h = h Q R^ s Y d ai W e 2i+ \ (3.3) 

where 7? > 0, \ 6\ > 1.154, the critical exponents for the 3-dimensional Ising model read j8 = 0.325, 
S = 4.80, a = 2- j3 (1 + 8) and a x = 1, a 3 = -0.76201, a 5 = 0.00804. The function g(6) follows 
from solving the differential equation 

2 

£ aa+i e 2 ' +1 (1-0 2 + 2j80 2 ) = 2(2 - a)dg(d) + (1 - d 2 )g'(6) (3.4) 

with integration constant g(d = 1) = 0.04242. Following these definitions, M(r,/i) shows the cor- 
rect critical behavior M(r = 0,h) ~ sgn(/i)|/i| 1 / 5 , M(r < 0,h = + ) ~ \r\& close to r = 0, h = 
when choosing Mo, /zo appropriately. In Fig. 3, we illustrate the mapping of the variables of the 




Figure 3: Mapping (r,h) (R, 9) of the parametric representation of G in |3(]] according to (3.2, 3.3). 
Left panel: solid curves denote lines of constant /; = 0.5 (—0.5) for the upper (lower) curve. Dashed curves 
depict lines of constant r > (r = 0.5, 1 , 2 from left to right), whereas dotted lines represent constant negative 
r = —0.5. Right panel: solid curves denote lines of constant R = 0.5, 1,2,3 from inner to outer ring, dashed 
curves lines of constant 0. For h < 0, is negative (9 = —1.1,-1 (here r = 0) and —0.5 from left to right); 
for h > 0, 9 is positive (9 = 1.1,1 (here r = 0) and 0.5 from left to right). For h = 0, 9 = ±1.154 for r < 
and 0=0 for r > 0. By construction, the parametrization is continuous for r > 0, whereas the discontinuity 
of the paramettization in the region of r < generates a first-order phase ttansition for Hb > when 
choosing the coordinate system (r,h) in the T - jis plane appropriately, cf. Fig. 4. 



3-dimensional Ising model (r,h) on the parametric representation (R, 9) of the Gibbs' free energy 
density G constructed in [28, 3C]. 

CP is located at r = and h = 0, cf. Fig. 4, and the coordinate system is orientated such 
that by construction the region of positive r denotes the cross over region, whereas the region 
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Figure 4: Visualization of the coordinate transformation between (T^/Ib) and {r,h). CP is located according 
to [||] on the estimated phase boundary T c (ij.b). The coordinate system (r 7 h) is located such that r = 0, h = 
at CP. The r-axis is defined to be tangential to T c (jXs) at CP such that r > denotes the cross over region, 
r < the region of first-order phase transitions. We put the /z-axis perpendicular to the r-axis. Note, that 
with increasing fig, r-axis and phase boundary increasingly deviate from each other. 



of r < defines the region of first-order phase transitions. The phase boundary is estimated by 



|, |3Tj] which agrees also fairly well 
175 MeV. 



T c (n B ) = T c (\ + jd(^r) 2 ^j with d = -0.122 according to 
with the curvature when solving for the flow equation of our QPM emanating at T c 

The singular part of the entropy density near the QCD critical point can be constructed from 
the dimensionless quantity 



Sc{T,pL B ) 



AT? + An 2 



B.c 



dG 
dT 



(3.5) 



Here, AT C and A/i B c approximate the extension of the critical region in T and jj, B directions, re- 
spectively, provided the CP is located at small jig, by quantifying the rapidity with which S c (T,jj, B ) 
changes across the r-axis that mimics the phase transition line. From Fig. 4 it also becomes clear 
that this phenomenological construction of CP features into the model is limited in accuracy by 
the increasing discrepancy between the r-axis and the actual phase boundary between hadronic and 
quark-gluon matter as }Xb increases. Thus, this procedure works best if CP is located in the flat 
region of the phase border line at small Hb- The difficult task is to map the QCD variables T and 



Hb on r and h and to evaluate (3.5) from (3.1). These relations were elaborated in the pioneering 



work [ 28 [ tne interested reader is referred to for the explicit form. 

In order to illustrate CP effects on the EoS, in particular on the pattern of isentropic trajec- 
tories, we consider here Nf = 2 and the CP located according to [||] at pL^ p = 360 MeV. Since 
we use r c (jU# = 0) = 175 MeV ([25], see also [32]) we find by solving the flow equation a 
slightly different T CP = 170 MeV (compared to []3|]). For the regular contribution s reg we employ 
( p| ) including 50,2,4,6(7") as calculated from our QPM with parameters as applied in [^3|, |2^] for 
describing c,(r) (cf. Fig. 2 for S2^{T))- When constructing s S i„ g from the dimensionless quantity 
S c (T,Hb) in (|3~3|), one important constraint to be satisfied is a vanishing net baryon density «g at 
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Hb = 0. From standard thermodynamic relations one finds 

n B (T,n B )-n B (0,n B )= C ^p^- dT' . (3.6) 
Jo a\l B 

Assuming a finite extension of the critical region around CP, the integration constant n B (0,n B ) is 
completely determined by the regular part of the thermodynamic potential. Obeying the correct 
dimension of an entropy density, we therefore choose s s i ng proportional to the second-order Taylor 



expansion coefficient c%{T) from ( |2.3[ ) (cf. Q230), 

s sing (T,IX B ) = ^Ac 2 (T)pL 2 B TtanhS c {T : pL B ). (3.7) 

Here, A denotes the relative strength between regular and singular contributions and tank S c binds 
S c (T,Hb) between —1 and +1 for convenience. 

In order to describe lattice QCD results, one has to adjust the remaining parameters AT C , A/j.b.c 
and A according to s/qcd = s reg + s s i ng , where for small jis lattice QCD results have already been 
successfully described by the regular part alone. A possible solution is choosing the parameters in 



(3.5, 3.7) as AT C = A/j.b. c = 2 MeV mimicking a small critical region and A = 0.1. This solution 
implies that bulk thermodynamic quantities like p(T,/j.b), e(T,H B ) or s(T,fi B ) are represented to a 
large extent by the regular contribution, i. e. by the QPM parametrization, as evident in Fig. 5 (left 
panel) for the isentropic trajectories [f33|]. Only for larger in the vicinity of the CP, the isentropic 
trajectories are slightly modified, however, leaving the general pattern unchanged. 

Increasing A or AT C and A^i B c and therefore the strength of the singular contribution or the ex- 
tension of the critical region, respectively, would change the isentropic trajectories also for smaller 



Hb- For instance, in [ ]33[ ] a simple toy model including schematically CP effects was discussed 
where for larger A > trajectories were bent into the opposite direction towards larger due to 
the presence of CP showing clear signatures of an attractor. Furthermore, trajectories on the right 



of CP exhibited the typical first-order behavior, cf. Fig. 5 (right panel). Note that in [ |28[ ] a model 
was constructed involving a hadronic low-temperature and a partonic high-temperature phase in 
which CP effects also change the pattern of the trajectories. These results significantly depend on 
the entropy density, in particular in the hadronic phase, but also on the strength parameter A and 



the extension of the critical region. However, lattice QCD results [ ]25fl seem to exclude such strong 
effects of CP. 

The small effect observed for small n B in Fig. 5 is understood by the numerical smallness 
of s S i n g (cf- Fig. 6 left panel) and its /^-dependence in comparison with s reg (cf. entropy density 
s(T,/i q = 0) in Figure 1 of [|20|]). At larger values of }Xb, the influence of CP increases, but with 
our focus on the small-/^ region we conclude that CP effects may be small for small net baryon 
densities, in particular on averaged hydrodynamics. Quark-meson model based studies combined 



with the proper-time renormalization group method [11] also point to CP effects concentrated on a 
fairly narrow region. 

In contrast, for close to CP, baryon number susceptibility %b = (dn B /dfi B )T or specific 
heat ch = T(ds/dT)^ B are strongly influenced by the QCD critical point. At CP and across the 
first-order phase transition line, c# diverges due to its singular part (cf. Fig. 6 right panel). In order 
to illustrate the effects on %b, we increase A = 1 which does not alter the qualitative behavior of 
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Figure 5: Influence of CP on the pattern of isentropic trajectories. Left panel: s/hb = 300, 45, 35 (from left 
to right) for the QPM without CP (solid curves, cf. Fig. 1) and with CP (dashed curves). For small /ig the 
influence of CP on the EoS is negligible. CP is located according to |^] on the estimated phase border line 
(dotted curve). Note that in our approach s reg and therefore s s i ng are based on a Taylor series expansion which 



is limited by its radius of convergence, /1b/T < 1 .8 [ 3 1 1 . It is therefore impractical to examine trajectories on 
the right of CP within this ansatz. In contrast, the full QPM is not hampered by these shortcomings [ [23| ] and 
it would, therefore, be conceivable to apply the full entropy density s = dp/dT from the QPM expression for 
p as regular part, defining s S i ng — As reg tanhS c - These considerations will be reported elsewhere. Here, we 
employ a toy model for illustration purposes in the right panel: Influence of CP on the isentropic trajectories 
s/ns — 50, 28 (from left to right) in a toy model. Here, s reg = 4c~oT 3 + ^c~2ligT, s s i„ g — ^^pigTAtanhSc 
with c = (32 + 2 17V/) n 2 / 1 80, c 2 = N f /2 and Nf = 2 where A — 0,0.5, 1 for dashed, thin and solid curves, 
respectively. For details cf. jp3[|. Dotted curve represents the r-axis. The isentropic trajectories display the 



typical pattern when crossing a first-order border line [ 34 1 signalling a finite latent heat. 



the observables. In Fig. ^ (left panel), we exhibit the baryon number density yib as a function of 
which shows the expected continuous behavior in the crossover regime and also for T > T c . At CP, 
however, the slope tremendously increases becoming infinite which turns into a finite discontinuity 
across the first-order phase transition line. Accordingly, its derivative %b as measure for fluctuations 
in jib diverges at CP and across the phase transition line whereas %b remains finite for temperatures 
T > T CP (cf. Fig. fright panel). Also, a suppression of baryon number fluctuations in the confined 
phase (solid, dashed and dash-dotted curves in Fig. |7| right panel) are evident in the vicinity of 
CP, whereas this is not the case in the deconfined phase (dotted curve). In contrast, the isovector 



susceptibility Xi is expected to remain finite when approaching CP []24]. Here, we examined the 



behavior of Xb for constant temperatures. In [ 35 ] it was argued that a verified signal of CP would 



require a non-monotonic behavior of Xb along the phase boundary T c {}JLb)- 

Other possibilities of including CP in singular contributions of thermodynamic bulk quantities 



are conceivable, e. g. s s i ng in (3.7) could be chosen proportional to rather than jU^T. Also, the 



inclusion of CP effects into different toy models were reported in [B3H. Different phenomenological 



approaches of constructing an EoS with CP behavior can be found in [36, 37]. However, a unique 



adjustment to lattice QCD data without better information or guidance seems not to be accessible. 
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Figure 6: Left panel: scaled singular part of the entropy density Ssing/T^ (3.7) as a function of T/T c for 
different pLg = 210, 360, 450 MeV (solid, dashed and dash-dotted curves, respectively). For small pig, s s t ng 
is negligible compared to s reg (cf. entropy density s(T,ji q = 0) in Figure 1 of [poj]). At CP the slope of 
s S i n g becomes infinite resulting in a diverging specific heat ch = T4f. For /ig > 360 MeV, s s i ng shows the 
significant discontinuity across the phase transition line of first-order. Right panel: corresponding singular 
part of the specific heat ch (same line code). 



4. Conclusions 



In summary we present here a phenomenologically guided construction to supplement a simple 
quasi-particle model by elements displaying the static criticality of the 3-dimensional Ising model, 
i. e. the universality class QCD belongs to. In doing so, we apply a procedure well probed in solid 
state and condensed matter physics, but which is restricted to a narrow region around the critical 
point, provided its existence. Having adjusted the pure quasi-particle model to available lattice 
QCD results, little space is left for pronounced effects related to the critical point. Of course, an 
divergence of susceptibilities at the critical point and the first-order phase transition beyond it (at 
larger chemical potentials) are found, but the modification of the equation of state remains rather 
modest, as evident in the isentropic curves in the T - Hb plane. It turns out that various ansatze 
are conceivable and that thermodynamic consistency restricts the parameter space. The latter fact 
leaves us with the impression that a readjustment of the quasi-particle model parameters, such that 
even for sizeable critical point effects the lattice QCD results are recovered, is hardly possible. 

One major prospect of observable consequences of the critical point is the occurrence of in- 
creased fluctuations. An adequate theoretical study has to rely on a more microscopical frame- 
work than offered by using the equation of state in standard hydrodynamics, say along the lines in 



36J, |38p. Nevertheless, the presented model may be used in further studies, e. g. considering hydro- 
dynamical expansion along neighboring paths crossing or not the first-order phase transition region 
or the critical point itself and analyzing sensible flow characteristics like elliptic flow and higher 
moments thereof. Another issue may be the study of critical point effects on penetrating probes, 
like photons and dileptons [39], i. e. the electromagnetic emissivity of hot and dense strongly 
interacting matter. 
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Figure 7: Left panel: baryon number density ng as a function of fig for different temperatures T = 
T CP - 1.25 MeV, T CP , T CP + 1.25 MeV and T = 176 MeV (dashed, solid, dash-dotted and dotted curves, 
respectively). For T > T CP , is continuously increasing, however showing an increasing influence of 

CP's presence by a tiny depression when approaching T CP ; at T = T CP , the slope becomes infinite at /i^ ', 
and for T < T CP , «b(Mb) exhibits the discontinuous behavior of a first-order phase transition when crossing 
the phase boundary. Right panel: corresponding baryon number susceptibility %b = as function of Hb 
(same line code). The divergence of %b{^b P ) * s clearly seen, while the spike for T — T CP — 1.25 MeV at 
/ig = 0.41 GeV is an artifact of the numerical differentiation at the density discontinuity. Note that due 
to the singular contribution we are restricted in our choice of parameters in (3.5, 3.7) in order to maintain 
thermodynamic self-consistency. The parameters here are chosen such that this possibility becomes visible 
in the region near Hg P ~, where Xb(^b) drops as a consequence of the depression in n^/ig). In principle, for 
an unfortunate parameter choice, the slope of «b(Mb) in the vicinity of CP can become negative resulting in 
a negative % B . 



Acknowledgment 

This work is supported by BMBF 06 DR 121, 06 DR 136, GSI, Helmholtz association VI and 
EU-I3HP. We thank M. Asakawa, F. Karsch, E. Laermann, A. Peshier, K. Redlich and J. Wambach 
for fruitful discussions. 



References 

[1] K. Rajagopal and F. Wilczek, At the frontier of particle physics, vol. 3, Ed. M. Shifman, World 
Scientific, p. 2061 [hep-ph/ 001 1333]. 

[2] D. H. Rischke, Prog. Part. Nucl. Phys. 52, 197 (2004) [nucl-th/0305030]. 

[3] Z. Fodor and S. D. Katz, JHEP 0203, 014 (2002) [hep-lat/ 01 6002]; JHEP0404, 050 (2004) 

[hep-lat/0402006]. 

[4] R. V. Gavai and S. Gupta, Phys. Rev. D 71, 114014 (2005) [hep-lat/ 04 12 035]. 

[5] M. A. Stephanov, Prog. Theor. Phys. Suppl. 153, 139 (2004); Int. J. Mod. Phys. A 20, 4387 (2005) 
and references therein [hep-ph / 04 02115]. 

[6] P. de Forcrand andO. Philipsen, Nucl. Phys. B 673, 170 (2003) [hep-lat /0307 02 0]. 



10 



Quasi-p article perspective on critical end-point 



Burkhard Kampfer 



[7] C. Schmidt, C. R. Allton, S. Ejiri, S. J. Hands, O. Kaczmarek, F. Karsch and E. Laermann, Nucl. 
Phys. Proc. Suppl. 119,517 (2003) [hep-lat/0209009]. 

[8] M. A. Halasz, A. D. Jackson, R. E. Shrock, M. A. Stephanov and J. J. M. Verbaarschot, Phys. Rev. D 
58, 096007 (1998) [hep-ph/9804290]. 

[9] P. de Forcrand and O. Philipsen [hep - 1 at / 6 7 1 7 ] . 

[10] O. Philipsen, Proc. Sci. LAT2005, 016 (2005) [hep-lat/051007 7]. 

[11] B.-J. Schaferand J. Wambach [hep-ph/0603256]. 

[12] Y. Hatta and T. Ikeda, Phys. Rev. D 67, 014028 (2003) [hep-ph/0210284]. 

[13] F. Karsch, E. Laermann and C. Schmidt, Phys. Lett. B 520,41 (2001) [hep-lat/0107020]. 

[14] R. D. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 (1984). 

[15] E. Laermann and O. Philipsen, Ann. Rev. Nucl. Part. Sci. 53, 163 (2003) [hep-ph/0303042]. 
[16] F. Karsch, Lect. Notes Phys. 583, 209 (2002) [hep-lat/010 6019]. 

[17] M. Asakawa, U. Heinz and B. Miiller, Phys. Rev. Lett. 85, 2072 (2000) and references therein 

[hep-ph/0003169]. 

[18] A. Peshier, B. Kampfer, O. P. Pavlenko and G. Soff, Phys. Lett. B 337, 235 (1994); Phys. Rev. D 54, 
2399 (1996). A. Peshier, B. Kampfer and G. Soff, Phys. Rev. C 61, 045203 (2000) 
[hep-ph/9 911474];Phys. Rev. D 66, 094003 (2002) [hep-ph/02 622 9]. 

[19] M. Bluhm, B. Kampfer and G. Soff, J. Phys. G 31, 1151 (2005) [hep-ph/0411319]. 

[20] M. Bluhm, B. Kampfer, R. Schulze and D. Seipt, Eur. Phys. J. C in print [hep-ph/ 60 8 053]. 

[21] M. Le Bellac, Thermal Field Theory, Cambridge University Press, Cambridge, 1996. 

[22] M. I. Gorenstein and S. N. Yang, Phys. Rev. D 52, 5206 (1995). 

[23] M. Bluhm, B. Kampfer and G. Soff, Phys. Lett. B 620, 131 (2005) [hep-ph/041110 6]. 

[24] C. R. Allton, S. Ejiri, S. J. Hands, O. Kaczmarek, F. Karsch, E. Laermann and C. Schmidt, Phys. Rev. 
D 68, 014507 (2003) [hep-lat/0305007]. C. R. Allton, M. Doring, S. Ejiri, S. J. Hands, O. 
Kaczmarek, F. Karsch, E. Laermann and K. Redlich, Phys. Rev. D 71, 054508 (2005) 

[hep-lat/0501030]. 

[25] S. Ejiri, F. Karsch, E. Laermann and C. Schmidt, Phys. Rev. D 73, 054506 (2006) 

[hep-lat/0512040]. 

[26] M. Bluhm, B. Kampfer, R. Schulze and D. Seipt [hep-ph/0608052]. 

[27] F. Karsch, E. Laermann and A. Peikert, Phys. Lett. B 478, 447 (2000) [hep- 1 at / 0002003]. 

[28] C. Nonaka and M. Asakawa, Phys. Rev. C 71, 044904(2005) [nucl-th/050 90 91]. 

[29] W. Gebhardt and U. Krey, Phasenubergdnge und kritische Phanomene, Friedr. Vieweg & Sohn, 
Braunschweig/Wiesbaden, 1980. 

[30] R. Guidaand J. Zinn-Justin, Nucl. Phys. B 489, 626 (1997) [hep-th/9610223]. 

[31] C. R. Allton, S. Ejiri, S. J. Hands, O. Kaczmarek, F. Karsch, E. Laermann, C. Schmidt and L. 
Scorzato, Phys. Rev. D 66, 074507 (2002) [hep- 1 at / 0204010]. 

[32] M. Cheng et al., Phys. Rev. D 74, 054507 (2006) [hep-lat/060801]. 



11 



Quasi-p article perspective on critical end-point 



Burkhard Kampfer 



[33] B. Kampfer, M. Bluhm, R. Schulze, D. Seipt and U. Heinz, Nucl. Phys. A 774, 757 (2006) 
[hep-ph/ 5 9 1 4 6]. M. Bluhm and B. Kampfer [hep-ph/ 511015]. 

[34] H. W. Barz, B. Kampfer, L. P. Csernai and B. Lukacs, Phys. Lett. B 143, 334 (1984). 

[35] K. Redlich, B. Friman and C. Sasaki [hep-ph/ 60 92 1 1]. 

[36] Y. Hama, R. P. G. Andrade, F. Grassi, O. Socolowski, Jr., T. Kodama, B. Tavares and S. S. Padula, 
Nucl. Phys. A 774, 169 (2006) [hep-ph/0510096]; AIP Conf. Proc. 828, 485 (2006) 
[hep-ph/0510101]. 

[37] C. Ratti, M. A. Thaler and W. Weise, Phys. Rev. D 73, 014019 (2006) [hep-ph/ 050 6234]; 
[nucl-th/0604025]. 

[38] K. Paech, H. Stocker and A. Dumitru, Phys. Rev. C 68, 044907 (2003 ) [nucl-th/0302013]. K. 
Paech and A. Dumitru, Phys. Lett. B 623, 200 (2005) [nucl-th/ 0504 03]. K. Paech et al., these 
proceedings. 

[39] K. Gallmeister, B. Kampfer and O. P. Pavlenko, Phys. Lett. B 473, 20 (2000) [hep-ph/ 9 9082 6 9]; 
Phys. Rev. C 62, 057901 (2000) [hep-ph/000 6134]. 



12 



